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SUMMARY 

A  numerical  technique  has  been  developed  to  compute  the  optimum 
spanwise  load  distribution  on  nonplanar  wings  of  arbitrary  shape. 
Complex  curved  wing  configurations  with  multiple  fences  can  be  analyzed 
with  this  technique.  M.  M.  Munk's  criterion  for  minimum  induced  drag 
was  used.  The  problem  is  solved  in  the  two-dimensional  Trefftz  plane. 
The  two-dimensional  shed  vortex  sheet  is  assumed  to  have  the  same  shape 
as  -he  nonplanar  wing,  from  which  it  has  been  shed. 

J.  L.  Lundry  of  the  McDonnell  Douglas  Corporation  found  an 'ingen¬ 
ious  solution  to  this  problem  by  computing  the  potential  flow  field, 
which  would  satisfy  Munk's  criterion.  Lundry's  method  requires  a 
Schvartz-Christoffel  conformal  transformation. 

The  method  developed  in  this  paper  is  different  in  that  the 
numerical  technique  does  not  require  a  conformal  transformation.  The 
vortex  sheep  in  the  Trefftz  plane  is  subdivided  into  2N  segments. 

Each  vortex  sheet  segment  is  assumed  to  have  a  linear  vorticity  distri¬ 
bution.  The  velocity  induced  at  N~Q  stations  is  determined  with  the 
Biot-Savart  law.  Because  of  symmetry  it  is  sufficient  to  compute  the 
velocities  in  one  half  of  the  vortex  sheet.  A  set  of  N-Q  equations 
has  been  derived  with  as  many  unknowns.  The  unknowns  are  the  strength 
of  the  shed  vortex  sheet  at  N-Q  stations.  Munk's  criterion  provides 
the  condition  matrix  for  the  magnitude  of  the  normal  component  of  the 
induced  velocity.  The  strength  of  the  shed  vortex  sheet  is  integrated 
to  obtain  the  optimum  spanwise  loading  on  the  nonplanar  wing. 
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The  technique  has  been  tested  on  nonplanar  wings  with  various 
dihedral  angles  and  locations  of  the  nonplanar  wing  sections.  The 
nonplanar  wing  results  are  presented  relative  to  those  of  a  planar 
wing  with  elliptical  loading.  Both  wings  elected  have  identical  lift 
and  total  wing  peripheral  length*  and  thus  equal  skin  friction  drag. 

It  is  shown  that  the  ratio  of  the  induced  drag  of  the  planar  wing  to 
that  of  the  nonplanar  wing  is  always  less  than  one.  The  results  are 
in  complete  agreement  with  those  obtained  by  Lundry.  However*  because 
Lundry  compared  wings,  with  equal  span  instead  of  equal  peripheral 
length,,  he  found  this  ratio  to  be  always  greater  than  1.0. 

In  conclusion  it  can  be  said  that  if  the  span  is  the  limiting  factor 
it  may  be  advantageous  to  use  nonplanar  wings  with  dihedral  angle  and 
fences.  However  if  the  wing  peripheral  length  is  limited,  then  the 
planar  wing  is  always  the  most  desirable  configuration ,  with  the  highest 
lift  over  drag  ratio. 
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abbreviations  as  defined  in  appendix 
induced  drag 

variable  along  vortex  sheet  segment 
induced  drag  efficiency  factor 
lift 

total  number  of  the  vortex  sheet  segments 
unit  normal  vector 

abbreviation  as  defined  in,  the  appendix 

total  number  of  ving  and  fence  tips  on  a 
semi span 

radius  vector 

radius  vector 

abbreviation  as  defined  in  the  appendix 
ratio  as  defined  in  Eq.  (27) 
dovnvash  velocity 

coordinates  of  the  vortex  sheet  segments 
unit  vector  in  downstream  direction 


Greek  Letters: 
'CH 

r 

r 

Y 


coefficient  defined  in  Eq,  (ll) 
coef f  'jient  defined  in  Eq.  (30) 
bound  vortex  strength 
•shed  vortex  sheet  strength 
angle  defined  in  Pig,  4 
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SYMBOLS  (cont'd) 


dihedral  angle 

non-dimenBional  spamd.se  coordinate 
air  density 


Superscripts  s 

'  indicates  the  side  of  the  vortex  sheet  at 

vhich  the  dcvmmsh  is  computed 

"  indicates  opposite  side  of  the  vortex 

sheet  from  vhere.  the  dcvmr&sh  is  computed 

*  properties  of  a  planar ,  elliptic.ally 

loaded  ving 


Subscripts : 

j  ■  indicates  anyone  of  the  Q  tip  sections 

on  the  semispan 

m  indicates  the  location  of  the  inducing 

vortex  sheet  segment 

n  indicates  the  location  at  vhich  the 

induced  velocity  is  computed 

o  indicates  the  vortex  sheet  centerline 

conditions 
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INTRODUCTION 

A  progressively  increasing  need  is  11600111108  apparent  for  aircraft 
with  STOL  capabilities.  The  forty  plus  versions  of  fixed  wing  air¬ 
craft  existing  todey,  with  both  VTOL  and  STOL  capabilities,  represent 
a  bare  initial  effort  to  satisfy  both  commercial  and  military  require¬ 
ments.  Commercial  nee's  include  downtovn-to-airport  transportation, 
while  those  of  the  military  are  for  operations  from  short  runways  such 
as  those  found  on  aircraft  carriers . 

STOL  capability  requires  unusual  wing  configurations,  often  using 
multiple  leading  edges  and  trailing  edge  flaps  with  boundary  layer 
control  to  obtain  high  lift  coefficients.  These  modifications  often 
require  many  fences  to  separate  wing  sections  with  different  lift 
coefficients.  The  frequent  use  of  horizontal  or  vertical  engine  nacelles 
adds  further  to  the  nonplanar  nature  of  the  wings.  Aircraft  on  Naval 
carriers  can  benefit  from  a  nonplanar  wing  configuration  so  as  to  redv.cc 
their  required  span  aud  permit  compact  storage. 

The  computer  technique  developed  in  this  paper  calculates  the 
optimum  loading  on  any  nonplanar  wing  of  given  complex  geometry  for  the 
minimum  lift  over  drag  ratio.  Such  a  computer  program  can  be  used  to 
determine  the  effect  of  size,  location,  and  dihedral  angle  of  the  fences 
employed  on  the  overall  lift  to  induced  drag  ratio  of  the  wing.  The 
computed  optimum  loading  can  determine  the  desired  wing  twist  end 
taper. 
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The  optimization  technique  is  based  on  Muuk's1  criterion  for 

,1 

minimum  induced  drag.  This  criterion  requires  that  the  cosine  of  the 
local  dihedral  angle  equal  the  normal  component  of  the  local  induced 
velocity  divided  by  a  constant. 


Figure  Is  Munkrs  criterion  f  r  minimum  induced  drag* 

For  all  planar  wings,  regardless,  of  taper,  aspect  ratio,  or  sweep, 
the  optimum  spanvise  loading  is  elliptical. 

Because  the  induced  drag  of  a  .ring  is  only  a  function  of  its 
spanvise  loading  and  not  of  the  aerjzct  ratio  it  is  desirable  to  derive 
the  optimum  spanvise  loading  in  the  Trefftz  plane,  which  is  located 
infinitely  far  downstream  of  the  wing,  where  all  chordwiBe  effects  can 
he  ignored.  The  two-dimensional  Trefftz  plane  is  perpendicular  to 
thb  .  ,>uble-inf inite  shed  vortex  sheet ,  which  is  assumed  to  remain 
undistorted  and  to  maintain  the  same  nonplanar  shape  as  the  wing. 

^•Superscript  numbers  denote  references. 
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Lundry2  has  presented  a  solution  to  ii i  problem  by  computing  the 
potential  flow  field  which  satisfies  Munk's  criterion.  By  superim¬ 
posing  a  uniform  crossflow  on  the  two-dimensional  vortex  sheet,  he 
cancelled  all  the  normal  components  of  the  downwash  velocity  and  thus 
satisfied  the  boundary  conditions  on  a  solid  surface.  Next  he 
computed  the  velocity  potential  around  a  solid  surface  in  the  shape 
of  the  nonplanar  wing  by  using  a  Schvartz-Christoffel  conformal  mapping 
technique.  Then  by  subtracting  the  velocity  potential  of  the  super¬ 
imposed  crossflow,  he  found  the  desired  velocity  potential  of  the 
induced  velocities.  The  strength  of  the  shed  vortex  sheet  was  found 
from  the  difference  in  the  upper  and  lower  tangential  velocity 
components.  Using  integration,  he  subsequently  computed  the  optimum 
spanwise  loading  on  the  wing. 


The  velocity  past  a  solid  surface  in  the  shape  of  the  nonplanar 
wing  is  zero  at  stagnation  poilnvs #  e.g. ,  the  root  of  a  fence  on  the 
wing,  and  is  infinite  for  the  flow  along  the  wing  where  there  is  a 
discontinuity  in  the  dihedral  angle,  such  discontinuities  will  he 
referred  to  as  bends.  Consequently  Lundry's  analysis  shews  that  the 
strength  y  of  the  vortex  sheet  segment  shed  by  the  fence  is  zero 
at  the  fence  root.  In  addition  a  vortex  sheet  segment  shed  by  either 
the  top  of  a  wing,  fence,  or  by  the  bend  in  the  wing,  has  infinite 
vortieity.  Such  segments  require  special  treatment  in  a  numerical 
calculation. 

Lundry's  analysis  covers  straight  line  segments  because  of  the 
nature  of  the  Schvartz-Christoffel  transformation,  but  others^  have 
used  mappings  which  have  transformed  curved  lines  to  straight  lines 
and  performed  a  similar  analysis.  In  references  2  and  3,  the  analysis 
wa3  performed  in  the  Trefftz  plane  and  an  undistorted  wake  assumption 
was  made.  A  recent  publication  by  Blackwell1*  describes  a  numerical 
technique  which  does  not  use  conformal  mapping.  His  technique  is  a 
special  form  of  lifting  surface  theory,  which  provides  information 
on  both  the  spanwise  and  chordwise  optimum  loading.  Experimental 
work  on  nonplanar  wings  has  been  reported  by  Roberts^. 


METHOD  OF  ANALYSIS 


The  vortex  sheet  shed  by  a  given  wing,  with  arbitrary  shape 
and  number  of  fences,  is  assumed  to  extend  to  infinity  without  any 
distortion.  Only  wings  with  an  optimum  loading  will  have  a  trailing 
vortex  sheet  which  satisfies  Munk's  criterion  at  any  spanwise 
location: 


CO.  ©„  =  - ZkA. .  (!) 

We=  constant  ' 

This  criterion  can  be  applied  to  the  spanwise  direction  at  the  wing 
or  at  any  other  normal  surface  downstream  of  the  wing.  If  the 
normal  surface  is  chosen  infinitely  far  downstream  of  the  wing  then 
one  deals  with  the  Treffiz  plane  and  the  problem  has  become  two- 
dimensional.  The  velocities  found  in  the  Trefftz  plane  are  induced 
by  a  two-dimensional  double-infinite  vortex  sheet  of  strength  y 
The  shape  of  the  shed  vortex  sheet  is  equal  to  that  of  the  nonplanar 
wing,  assuming  an  undistorted  wake. 

Identifying  the  semi  wing  total  peripheral  length  by  y  max  *  1.0 
and  the  bound  vorticity  of  the  wing  by  P  a=  P  (y')  ,  then  the 

strength  of  the  shed  vortex  sheet  is: 


If  the  seai-«wfee)s  she«t  total  periphtuv'.  <.v .  w»  •  lsax. 


"  :  J 


iatJ  N 


of  e^ual  length  A-v?  as  aho wn  in  Fig*  3,  then; 


A'*?  » 


N 


(3) 


Define  the  number  of  wing  and  fence  -tips  on  the  semispan, 
if  each  tip  has  N(j)  vortex  shoot  segments s  where  1,< 
then  the  total  number  of  vortex  sheet  segments  is: 


it 


4. 


N 
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N(j) 


An  example  it  shown  here  with  Q  »  5  and  N  **  26. 
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Figure  3:  Example  of  an  arbitrary  nonplanar  wing  shape 
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Hunk's  criterion  will  "be  satisfied  at  N-Q  stations  n^j, 
where  n  varies  from  1  through  1.  The  n^j  stations  are  located 

on  the  right  hand  side  of  the  vortex  sheet;  hecauae  Munk's  criterion 
is  automatically  satisfied  on  the  left  hand  side  owing  to  symmetry. 

The  normal  component  of  the  velocity  induced  at  j  by  the 
entire  shed  double  Infinite  vortex  sheet  is  V/j_  .  This 

velocity  component  is  composed  of  the  sum 


Ji&> 


W 


Nrfc-1 


A  W 


i-  n(i),  xa(j) 


(0 


where  AW,  ,  is  the  velocity  induced  at  n/ <  \  by  both 

JLn(3),m(i)  * 

the  left  hand  and  the  right  hand  vortex  sheet  segments  located 
between  stations  my )  and  m(  .  The  strength  of  these  vortex 
sheet  segments  is:  y  ("*})  for  the  right  hand  side  and  -  y(^) 

for  the  symmetrically  opposite  segment  on  the  left  hand  aie.  Using 
the  Biot-Savart  law  and  leaving  off  the  subscript  ( j )  gives 


■n,w 
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Where  r  and  R  can  be  expressed  as  functions  of  the  known  coordin¬ 
ates  of  the  vortex  sheet  segments: 
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R*,,m  =  (**  +  ^  eos  *0  C<5S  0I»  +  h  SIN  "  &n)  a,N  (f  *) 

-y^siw  (at) 


R  is  positive  for  lj/  <,  J  and  negative  for  y  |  -2L 
while  r2  is  positive  and  can  be  written  as: 


* K  &*-x»j *(tn  *  *>  5,“  (9  «) 

(*£,  *)  =  C" 11  ®».  -  *»)  +0™+  K  e»  -  %)  (9  l>) 


To  be  able  to  integrate  Eq,.  {7}  one  has  to  assume  a  vortex  distri¬ 
bution  along  the  segment.  Good  accuracy  has  been  obtained  by  assuming 
a  linear  vortex  distribution,  so  that  y  becomes  a  function  of  the 
two  unknown  values  at  the  ends  of  each  segment  Y  and  Y  : 

»**  l-VH*! 


r-  r» 


Jh—  r  y 
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Ccsabining  Eqs.  (7)  and  (10)  gives 
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Equations  (12-15)  are  integrated  as  shown  in  appendix  A  to  give 
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The  symbols  A,  B,  C,  etc.,  are  functions  of  the  Known  coordinates 
and  dihedral  angles  of  the  vortex  sheet  segments  (see  Eqs.  (38-1*7)  in 
Appendix  A).  Equations  (16-19)  are  the  general  forms  of  the  coefficients 
in  Eq.  (ll)  except  for  the  special  case  that  m  =>  n-1  or  m  »  n.  For 
m  =  n  one  finds  <*.  =.  00  and  for  m  «  n-1  one  finds  <*,  = 

This  is  easily  understood  when  it  is  realized  that  under  these  special 
conditions  the  contribution  of  the  right  hand  side  of  the  vortex  sheet 


computed  at  the  end  of  vortex  sheet  segments  of 

-  tn  =•>) 

finite  strength  and  and  then  the  velocity  becomes 

infinite.  This  problem  can  be  avoided  when  Lw,  1 

V'tM'flL 

*-  *  J  °<i  >  flfi 

is  computed  for  both  adjacent  vortex  sheet  elements  simultaneously  as 
shown  in  Fig.  ll*  (see  Appendix  B).  The  results  of  such^ special  treatment 
are; 

f 


for  m  «  n-1 


C*  SS  —  _CO*  £  Qw-gQ 

/  2L7T 


=sz- 
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for  m  «  n 


cos  A. 
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,  9(  JSS  4- 
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These  equations  have  been  derived  in  appendix  B  (see  Eqs.  (53-56)). 


12 


Another  difficulty  arising  in  the  application  of  Eq.  (11)  occurs 
when  m  +  1  «  or  Ym+i  s  which  is  the  vorticity  shed 

by  either  a  wing  tip  or  a  fence  tip.  This  vorticity  is,  of  coura?. 


infinite  and  when  used  in  Eq,  (li)  would  make  AW x  °° 

This  problem  has  been  overcome  by  assuming  that  y  shed  by  the  non- 
planar  wing  goes  to  infinity  at  at  the  same  rate  as  the  vorticity 
y  shed  by  an  elliptically  loaded  planar  wing  becomes  infinite  at 
the  edge.  Defining  the  semi span  of  the  planar  wing  by  ®  1.0  and 

the  strength  of  the  bound  vortex  at  the  centerline  by  H,  gives: 


_ .  * 

r  i 
°  0-?T 


and  y  at  the  distance  Z  Av  from  the  edge  is  given  by: 
'N-a  ' 


:  -r^.. 

'-2  0 -?/.,)* 

*  V*  ( 1  - 


r*=r*  ( 1  - 

°  v. 


0-71) 


The  normal  component  of  the  downwash  velocity  induced  at  station  R-l, 
adjacent  to  the  tip  is  most  sensitive  to  the  accuracy  of  the  vorticity 
at  the  tip.  If  the  downwash  velocity  at  N-l,  induced  by  the  vortex 
sheet  segment  between  N-l  and  N  is  computed*  the  result  will  be  infinite. 


is  evaluated  at  the  end  of  a  vortex  sheet 


This  is  because  AW,  „  is  evaluated  at  the  end  of  a  vortex  sheet 

.  r  ,  , 

with  finite  strength.  However,  if  A  W,  is  induced  by  both 

L  ■*-#-»  L 

adjacent  vejtex  sheet  segments,  then  this  normal  velocity  vill  be  finite, 


1  *■*  '  |  * 

with  finite  strength.  However,  if  A  W, 


induced  by  both 


as  shown  in  Eq.  (52)  (see  Appendix  B). 


13 


Ono  can  compute  [A  ^  induced  by  the  two  adjacent 

segments  of  strength  y*  as  shed'by  an  ulliptically  loaded’  planar 
ving  by: 


replace  y  by  cos  z  and  dy  by  -  sin  z  dz 


Compare  j^AW^ J  ^  to  the  expression  found  for  the  velocity 

induced  by  two  adjacent  vortex  sheet  segments  with  a  linear  vorticity 

distribution  as  was  derived  in  Eq.  {52)  (see  appendix  B)  and  which  is 
repeated  here  as: 


AW,  —  Jn(3)~s.  /  Yn(J)  cOS  '/(  6  —  0  (oe,\ 

iN(5j-l  ~  2  TT  V  Y  7  \26) 

U„c<}  lVy**»-*  '  J 

If  the  Induced  velocity  6W*  ]  is  set  equal  to  Lv/, 

eta  ,*  1  J<x,  L  -‘•NCi)-U«,,‘<j 

end  the  vorticity  y^  *qual  YN(i,-2  then  an  equivalent 
finite  value  for  the  tip  vorticity  yr  is  found  as  a  function  of 

the  vorticity  . 


Equating  Eq.  (25)  to  Eq.  (26)  gives  the  equivalent  finite  tip 
vorticity  yN  (i)  as: 

Tn  (3)  —  i 


where  T  is  given  by: 


Equation  (28)  is  quite  insensitive  to  the  size  of  the  increment  used. 
For  example: 


-1» 

*7  =  10 

gives 

T  =  5.986 

-3 

Ay  =  io 

gives 

T  «  5.991 

-2 

Ay  =10 

gives 

T  »  6.0M+ 

-l 

-dy  =10 

gives 

T  »  6.  Tilt 

7h-i  "  7n  -  ^7 
7n-£  ®  >  -  a  *7 


elliptically  loaded  wing 
shed  vorticity  y“* 


shed  vortex  sheet 


Figure  5:  Graphical  presentation  of  the  finite  tip  strength  approximation. 


All  the  terms  in  Eq.  (11)  have  now  been  defined  and  the  next  step 
is  to  solve  the  optimum  vorticity  distribution  in  the  shed  vortex  sheet. 
Applying  Munk's  criterion  from  Eq.  (1)  gives: 
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Using  the  expression  for  the  induced  velocity  from  Eq.  (11)  gives: 
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This  equation  can  be  rewritten  in  the  form: 
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where  the  coefficients  A  are  as  given  below,  with  subscript 
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Equation  (30)  represents  a  set  of  H-Q  equations  for  N-Q  unknown  values 
of  .  The  solution  of  this  set  of  linear  algebraic  equa¬ 

tions  is  obtained  with  the  use  of  the  digital  computer  subroutine 
DECOMP, ^  The  values  of  1  thus  obtained,  represents  the 

dimensionless  optimum  strength  distribution  in  the  shed  vortex  sheet. 

To  compute  the  corresponding  spanvise  distribution  of  the  bound 
vorticity  P  one  has  to  integrate  Eq.  (2)  numerically.  First, 
the  integration  is  performed  along  the  fences  starting  with  tip  N^) 


wing  root 


This  integration  is  carried  out  along  the  wing,  starting  at  the 
wing  tip  and  proceeding  to  the  centerline.  Every  time  a  fence 
root  is  encountered  the  bound  vortex  or.  the  wing  experiences  a  step¬ 
like  increase  of  magnitude  P($) 


wing  root 


The  location  of  position  n  is  shown  in  Fig.  6 
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RESULTS 


The  optimum  loading  computed  with  this  method  can  he  uacd  to 
calculate  the  corresponding  lift  and  induced  drag  or.  the  wing.  The  lift 
L  is  proportional  to  the  vertical  component  of  the  force  obtained 
from  the  vector  product  between  V<*,  and  P  and  is: 

W  *.('-0  c  os  9  cj  (») 

The  Induced  drag  3);  is  proportional  to  the  vector  product  w  x  P 

"jzi  i  ^ 

using  w  x  I  —  x  Tj-  W0  c«>s  qQy\  x  P} 

3>j  S  j  Poo  rjcos  Qdy  (34) 

"7im  x 

The  lift  over  induced  drag  ratio  for  wings  with  optimum  loading  is  then: 


L  _  V, 


To  compare  the  ratio  between  an  elliptically  loaded  planar 

wing  and  an  optimum  loaded  nonplanar  wing,  one  should  specify  that  both 
wings  must  have  equal  lift,  and  use  Eq.  (35)  to  define  an  induced  drag 
efficiency  factor  k  as: 
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planar 
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nonplanar 
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The  following  argument  shows  why  the  efficiency  factor  k  is  less  than 
unity  when  the  two  wings  being  compared  have  equal  total  peripheral 
length.  Eoth  wing3  are  chosen  to  have  equal  lift  or: 
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This  inequality  is  given  by  the  mean  value  theorem. 

fT 

Defining  an  average  bound  vorticity  I  = - .  j  [~i( 
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r  planar  ^  nonplanar 

If  one  defines  an  average  shed  vorticity  Y  ss  ~ —  J  Y  ^7  then 
according  to  Eq.  (2)  I  'v  y  •  .  From  Eq.  (30 )  it  follows 

that  Y  wd  or  combining  these  results  We.  %v»ax  • 

Assuming  the  proportionality*  constant  to  be  of  the  same  order  of 
magnitude  for  both  the  planar  and  the  nonplanar  wing  one  obtains: 
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This  is  the  definition  of  the  efficiency  factor  used 
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The  effect  of  a  dihedral  angle  on  the  efficiency  factor  k  ^  1.0 
is  shown  in  Pig.  7.  Both  planar  and  nonplanar  wings  have  the  same 
total  peripheral  length  and  consequently  the  same  area  and  skin 
friction  drag.  The  nonplanar  wing  has  a  shorter  span  hut  its  perfor¬ 
mance  is  worse  because  of  the  lower  lift  over  induced  drag  ratio. 

Seemingly  contradicting  results  are  found  in  Ref.  2,  where  the  effi¬ 
ciency  factor  k  1.0,  because  Lundry  has  compared  the  nonplanar  wing 
with  a  planar  wing  of  equal  span  instead  of  equal  peripheral  length. 

Under  these  conditions  the  area  and  the  skin  friction  of  the  nonplanar 
wing  is  larger  than  that  of  the  planar  wing  but  its  performance  is  better 
because  it  has  a  higher  lift  over  induced  drag  ratio  than  the  planar  wing. 
The  optimum  loading  computed  in  this  paper  is  in  exact  agreement  with 
the  results  obtained  by  Lundry  in  Ref.  2,  except  the  efficiency  factor 
k  is  defined  differently  (see  Pig.  7). 

In  conclusion  it  can  be  said  that  if  the  span  is  the  limiting  factor 
then  it  may  be  advantageous  to  use  nonplanar  wings  with  dihedral  angle 
and  fences.  However,  if  the  total  peripheral  wing  length  is  limited, 
then  the  planar  wing  is  always  the  most  desirable  configuration  with  the 
highest  lift  over  drag  ratio. 

The  numerical  solution  was  tested  for  a  simple  nonplanar  wing  with 
a  planar  center  section  and  a  dihedral  angle  ©  ,  which  is  constant 

for  the  outer  portion  starting  at  y  Q  .  The  optimum  vorticity 
distribution  in  the  shed  vortex  sheet  is  shuwn  in  dimensionless  form 
in  Figs.  8  and  9.  There  is  a  sharp  increase  in  y  near  the  bend  in 
the  wing,  H  »  100  in  this  calculation.  The  ratio  between  the  bound 
vcrticity  P  and  its  average  value  F  is  shown  in  Figs.  10  and  11. 

The  effect  of  the  dihedral  angle  on  the  required  increase  in  the  optimum 
loading  is  shown  in  Figs.  12  and  13. 


Figure  8 :  Optimum  vor 
in  the  Ireff ft 


Figure  9 :  Optimum  vort 
in  thelrefftz 


* 


Figure  1 2 :  Relative  bound  vorticity 
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CONCLUSIONS 

The  results  of  this  numerical  technique  have  been  compared  to  the 
conformal  mapping  solution  as  presented  hy  Lundry2  in  his  configuration 
5  and  excellent  agreement  was  obtained. 

The  numerical  solution  used  in  this  program  takes  little  computer 
time  and  is  suitable  for  computer  optimization  of  fence  position,  size, 
and  dihedral  angle.  Only  the  {3  matrix  coefficients  that  belong  to 
the  fence  in  question  need  to  he  recomputed  for  each  configuration 
change. 

The  step  size  Ay  need  not  be  constant  *  ^  as  was  used 
here.  Wherever  appears  in  the  equations  it  can  he  replaced  by 

the  local  magnitude  of  Ay  used. 

The  input  data  for  the  computer  should  include  the  total  number 
of  vortex  sheet  segments  N  and  their  coordinates  X^  ,  and 

dihedral  angle  9^  .  In  addition  A  y^  must  be  specified  for 

each  station  n,  if  the  segments  are  not  of  equal  length. 

It  is  noteworthy  to  realize  that  Lundry2  compared  his  nonplanar 
wing  to  an  elliptically  loaded  planar  wing  with  the  same  span  and 
found  the  induced  drag  efficiency  factor  k  }  1.0.  Consequently,  when 
the  span  is  the  limiting  factor  it  is  advantageous  to  use  a  nonplanar 
wing.  In  this  analysis,  the  nonplanar  wing  is  compared  to  an  ellipti¬ 
cally  loaded  planar  wing  with  the  same  peripheral  length.  The  induced 
drag  coefficient  k  was  found  to  be  ^  1.0,  and  consequently  if  the 
total  peripheral  length  is  a  limiting  factor  it  is  advantageous  to  use 
a  planar  wing.  Either  method  can  be  used  to  find  the  optimum  loading 
of  a  nonplanar  wing. 
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APPENDIX  A 


The  integration  of  the  coefficients 
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Equations  8a  and  8b  ,  can  be  -written  as: 
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APPENDIX  A 

Equations  9a  and  9b  can  be  rewritten  as: 


where  =  a  coj  +  2  sin 

and  +  (j^-  fcj)* 

where  =  2.  +  +2  Sin  jj^) 

and  ^-»v»ii  =  (**»* 

Inserting  these  expressions  in  Eqs.  12-15  gives: 


(?2) 

(43) 
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APPENDIX  A 

2 

The  denominators  are  of  the  type  ah  +  bh  +  c  =  H  ,vith  4ac  ^  b 
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APPENDIX  B 

The  special  case  is  given  where  m  »  n  and  m  »  n-1.  Under  this 
special  condition  the  normal  component  of  the  induced  velocity  at  n 
is  evaluated  at  the  end  of  a  vortex  sheet  segment  with  finite  strength 
VV  ;  consequently  the  contribution  to  iW,  from  the 

I™  A 

adjacent  segments  on  the  right  hand  side  of  the  vortex  sheet,  as  given 
by  the  coefficients  o(t  and  g(z  will  be  positive  and  infinite  for 
man  where  o( .  =s  +©0  and  be  negative  and  infinite  for  m  =  n-1 

where  =.—  00  This  infinite  velocity  can  be  avoided  if 

j»  «y\— f 

the  integration  is  carried  out  over  both  vortex  sheet  segments  simul¬ 
taneously  as  shown  in  Fig.  14. 


Figure  llf?  Vortex  sheet  segments  for  m  =  n  and  m  =  n-1. 


When  the  dihedral  angle  ^  0^,  ,  then  there  is  a  bend  in  the 

vortex  sheet  and  the  corresponding  exact  solution  of  ^8n(i~  00  • 
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In  order  to  smooth  out  the  tends  and  avoid  the  necessity  for  ^  to 
go  to  infinity,  one  can  choose  the  unit  normal  vector  along  the 
average  normal  to  the  two  segments  as  is  shown  in  Fig.  ll>. 

Applying  the  Biot  Savart  lav  as  in  Eq.  (8)  gives  the  velocity 

♦ 

component  induced  at  n  by  both  neighboring 

L  J  "*1  J  W|  j 

vortex  sheet  segments,  hut  not  including  the  contribution  from  the 
left  hand  side  of  the  vortex  sheet  as  given  by  0<x  and  &(s 


Linearizing  the  vorticity  distribution  as  in  Eq.  (10 )  produced  a  dis¬ 
continuity  in  y  at  station  n.  This  vorticity  distribution  can  be 
replaced  by  one  continuous  function  over  both  segments  y,  »  hy 
one  distribution  ya  which  goes  to  zero  at  station  n  so  that 
r*ri  +  n  as  shown  in  Fig.  lb. 
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As  a  result  of  this,  one  can  vrite  vith  the  use  of  Eq..  (11),  the  con¬ 
tribution  of  the  right  hand  side  of  the  vortex  sheet  as: 
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